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Dear Editor, 



This is the cover letter for the resubmission of our manuscript named 'Pair 
production of neutralinos via gluon-gluon collisions '(coded DN7258), writ- 
ten by Jiang Yi, Ma Wen-Gan, Han Liang, Yu Zeng-Hui and H. Pietschmann, 
to Physical Review Dl as regular article. 

^ . We did some modifications in our manuscript according to the referee's 

suggestions, 1. For the language problems: (1) The misspelled words are 
changed. (2) The symbol '+ X' is added for the gg parton states in the 
subtitle of section 2. (3) The bar over the sum symbol in formula (2.9) is 
changed, and the symbol for the luminosity in formula (2.10) and (2.11) is 
adjusted. 

2. For the more serious shortcomings: (1) The typing error in neutralino 
• mass matrix (2.2) is corrected. (2) The argument about divergences on 

page 5 is omitted. The word 'evolution' is substituted by the word 'calcu- 
lation'. (3) The illustration of 'the exchange of incoming gluons diragrams' 
is added. (4) The argument using Furry's theorem is changed. (5) The dis- 
cussion on the Higgsino-like neutralinos pair production is added in the text 
and a relevant figure named 'Fig. (3)' is also added. (6) The references of 
renormalization group routine we used to calculate our low energy scenario 
from mSUGRA and the width of the intermediate Higgs bosons are added 
in section 3. (7) Some discussions of the curves about the cross section of 



neutralino pair production as the function of tan (3 are added. 

3. For the additional questions: We use the newer sets of MRST of gluon 
densities to re-calculate the cross section and we found the newer results 
are only very little different between the MRST (mode 2) scheme and the 
older MRS (set G) scheme. This is addressed at the end of section 3. We 
also add the discussion of the uncertainty of scale variation in the text and 
add a figure named 'Fig. (2)'. 

We are very grateful for the referee's useful suggestions. 
With our best regards. 



Sincerely yours, 



Ma Wen-Gan 
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Abstract 

The production of a neutralino pair via gluon-gluon fusion is studied in 
the minimal supersymmetric model(MSSM) at proton-proton colliders. 
The numerical analysis of their production rates are carried out in the 
mSUGRA scenario. The results show that this cross section may reach 
about 80 femto barn for X1X2 P a i r production and 23 femto barn for 
X2X2 P & i r production with suitable input parameters at the future LHC 
collider. It shows that this loop mediated process can be competitive 
with the quark-antiquark annihilation process at the LHC. 
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1. Introduction 

The Standard Model(SM) is a successful theory of strong and electroweak interac- 
tions up to the present accessible energies. The hierarchy problem suggests that in principle 
the SM is the low energy effective theory of a more fundamental one. Supersymmetry (SUSY) 
is presently the most popular attempt to solve the hierarchy problem of the SM, where the 
cancellation of quadratic divergences is guaranteed and hence any mass scale is stable under 
radiative corrections. The most favorable candidate for a realistic extension of the SM is the 
minimal supersymmetric standard model(MSSM). In the MSSM, a discrete symmetry called 
R-parity|| is kept in order to assure baryon and lepton number conservations, because the 
gauge-coupling unification supports conservation of R-parity@]. It implies that there exists 
an absolutely stable lightest supersymmetric particle(LSP). Proper electroweak symmetry 
breaking induces the right properties of the LSP being a natural weakly interacting cold 
dark matter candidate, which can explain many astrophysical observations||. In most cases 
the LSP in the MSSM may be the lightest Majorana fermionic neutralino. The neutralinos 
are mass eigenstates which are model dependent linear combinations of neutral gauginos 
and higgsinos0. They are determined by diagonalizing the corresponding mass matrix. 
In the MSSM the mass matrix depends on four unknown parameters, namely /i, M 2 , Mi, 
and tan/3 = t^/fi, the ratio of the vacuum expectation values of the two Higgs fields, fi 
is the supersymmetric Higgs-boson-mass parameter and M 2 and M\ are the gaugino mass 
parameters associated with the SU{2) and U{\) subgroups, respectively. 

Direct search of supersymmetric particles in experiment is one of the promising tasks for 
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present and future colliders. The multi-TeV Large Hadron Collider (LHC) at CERN and the 
possible future Next Linear Collider(NLC) are elaborately designed in order to study the 
symmetry-breaking mechanism and new physics beyond the SM. When the LEP2 running 
will be terminated, the hadron colliders Tevatron and LHC will be the machines left in 
searching for supersymmetric particles. Therefore it is necessary to give a proper and full 
understanding of the production mechanisms of supersymmetric particles at hadron collid- 
ers. If supersymmetry really exists at TeV scale, SUSY particles should be discovered and 
it will be possible to make accurate measurements to determine their masses and other pa- 
rameters of the Lagrangian at the LHC, and then we will have a better understanding of the 
supersymmetric model]?]]. We know that there are several mechanisms inducing the produc- 
tion of a chargino/neutralino pair at hadron colliders. One is through the quark- ant iquark 
annihilation called Drell-Yan process, and another is via gluon-gluon fusion. Although the 
neutralino pair production via gluon-gluon fusion is a process induced by one-loop Feynman 
diagrams, the production rate can be still significant due to the large gluon luminosity in 
hadron colliders. The direct production channels of chargino/neutralinos and sleptons at 
the hadron colliders Tevatron and LHC pp/pp — ► XiXj + X and IV + X via quark- ant iquark 



annihilation are investigated by Beenakker et al||. A recent work [0 showed that the 
one-loop process of the lightest chargino pair production via gluon-gluon fusion at the LHC 
can be considered as part of the NLO QCD correction to the quark-antiquark annihilation 
process. The neutralino pair production via quark-antiquark annihilation at the LHC was 
also considered by Han et al.lj. The measurement of CP violating supersymmetric phases 
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in the chargino and neutralino pair production at the NLC was investigated by Barger et 



al||10||. The lower mass limit of 29.1GeV on the lightest neutralino is obtained experimen- 
tally by analyzing the DELPHI results with the assumption that M 1 /M 2 ^0.5,tan/3 = 1, 
fi = -62.3GeV and M 2 = AG.OGeV^. 

In this paper we concentrate on the direct neutralino pair production via gluon-gluon 
collisions at the LHC in the framework of the MSSM with complete one-loop Feynman 
diagrams. The numerical calculation will be illustrated in the CP conserving mSUGRA 
scenario with five input parameters, namely mu2, uio, Aq, fi and tan/3, where TO1/2, mo an d 
A are the universal gaugino mass, scalar mass at GUT scale and the trilinear soft breaking 
parameter in the superpotential respectively. From these five parameters, all the masses 
and couplings of the model are determined by the evolution from the GUT scale down to 



the low electroweak scale [|13|]. The paper is organized as follows: In section 2, we introduce 
the relevant features of the model and the analysis of the cross section in this work. In 
section 3, we discuss the numerical results of the cross sections and finally, a short summary 
is presented. 

2. The Calculation of pp -> gg + X -> XiX°j + x 

In the MSSM the physical neutralino mass eigenstates Xi = 1)2,3,4) are the combi- 
nations of the neutral gauginos, B and W 3 , and the neutral higgsino, H® and H®. In the 
two-component fermion fields ip® = — 2A 3 , ipjjo, ^i/o)0, where A' is the bino and A 3 is 

the neutral wino, the neutralino mass term in the Lagrangian is given by 
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where the matrix Y reads 



(2.1) 



Y = 



( Mi — m z sin 8w cos (3 m z sin Q w sin /3 \ 

M 2 cos 6w cos /3 — cos 6w sin /3 

— mz sin 6^ cos /3 mz cos 6^ cos /3 — /i 

\ mz sin9w sin/5 — cos 9w sin /3 — /i / 



(2.2) 



Mi, M 2 and /i can be complex and introduce CP- violating phases. By reparametrization 
of the fields, M 2 can be real and positive without loss of generality. The matrix Y is 
symmetric and can be diagonalized by one unitary matrix iV such that Np = N*YN + = 
diag(m^o,m^o,m^o,m^o) with the order of m^o < m^o < m^o < m^. Then the two- 
component mass eigenstates can be 



X-=N^l i,j 1 1. 



(2.3) 



The proper four-component mass eigenstates are the neutralinos which are defined in terms 
of two-component fields as 

(2.4) 



a ) (<=i,...,4), 



and the mass term becomes 



Xi 



(2.5) 



where Mj are the diagonal elements of N D . 

The neutralino pair via gluon-gluon collisions can only be produced through one-loop 
diagrams in the lowest order. The calculation for this process can be simply carried out 



by summing all unrenormalized reducible and irreducible diagrams and the result will be 
finite and gauge invariant. In this work, we perform the calculation in the 't Hooft-Feynman 
gauge. The generic Feynman diagrams contributing to the subprocess gg — > i n t ne 

MSSM are depicted in Fig.l, where the exchange of incoming gluons in Fig.l(a.l ~ 3) 
and Fig.l(c.l ~ 2) are not shown. All the one-loop diagrams can be divided into three 
groups: (1) box diagrams shown in Fig.l(a.l ~ 3), (2) quartic interaction diagrams in 
Fig.l(b.l ~ 2), (3) triangle diagrams shown in Fig.l(c.l ~ 2). The Z° boson interme- 
diated s-channel diagrams with quark or squark loops, as shown in Fig.l(b.2), Fig.l(c.l) 
and Fig.l(c.2), cannot contribute to the cross section. That can be explained in two fields. 
Firstly, the vector component of the Z° boson wave function does not couple to the initial 
gg state as the result of the Laudau-Yang theorem. Secondly, the CP-odd scalar component 
of the Z° boson does not couple to the invariant CP-even XiXj state. We should mention 
that there are also some diagrams not contributing to the process, which we do not draw 
them in Fig.l (for example the s-channel diagrams with trilinear gluon interactions). Since 
the vertices of A°(G°) — q — q vanish |14|, there is no diagrams with a triangle squark loop 
coupling with an A or G° Higgs boson. 

We denote the reaction of neutralino pair production via gluon-gluon collision as: 

g(Ps, fj)g(pi, v) — ► x° (Pi)x?(P2)- (2-6) 
The corresponding matrix element for each of the diagrams can be written as 

M = M S + M i + M a (2.7) 



where Ai s , hA l and Ai u are the matrix elements from s-channel, t-channel and u-channel 
diagrams, respectively. The amplitude parts from the u-channel box and triangle vertex 
interaction diagrams can be obtained from the t-channel's by the following: 

M u = ^_iy ijM i^^^ p3 ^ p4 ^^ u ^ (2.8) 

where the indices i and j are for the final neutralino. When i = j, there is a relative minus 
sign due to the Fermi statistics, which requires the amplitude to be antisymmetric under 
interchange of the two final identical fermions. The cross section for this subprocess at 
one-loop order in unpolarized gluon collisions can be obtained by 

Hs,g 9 - m = leb^)* C di E | - M|2? (2 - 9) 



where = 1/2 (m|. + rri-,. — s) ± y^m?. + m?. — s) 2 — Arn^.rn^. . The factor (\) Sij is due 
to the two identical particles in the final states. The bar over the sum means average over 
initial spins and colors. 

With the results from Eq.(2.9), we can easily obtain the total cross section at pp collider 
by folding the cross section of the subprocess a(gg — > XiX'j) with the gluon luminosity, 

r 1 dC 

<T( PP ^gg + X^x°x j +X)= dr—^a(gg - x° X° ,at S = rs). (2.10) 

^(m-o+m- ) 2 /s ar 



where y/s and denote the proton-proton and gluon-gluon c.m.s. energies respectively 



and is the gluon luminosity, which is defined as 

dC„ Q f 1 dx\ 



(2.11) 



Here we used r = X\X^ the definitions of X\ and X2 are from Ref.p|. In our numerical calcu- 
lation we adopt the MRST(mode 2) parton distribution function f g (xi,Q 2 ) [TO, and ignore 
the supersymmetric QCD corrections to the parton distribution functions for simplicity. The 
factorization scale Q was chosen as the average of the final particle masses 1(^1% + m Xj)- 
The numerical calculation is carried out for the LHC at the energy 14TeV. 

3. Numerical results and discussions 

In this section, we present some numerical results of the total cross section from the full 
one- loop diagrams involving virtual (s) quarks for the process pp — » gg + X — > XiXj + X, 
respectively. The input parameters are chosen as: m t = 173.8 GeV, mz = 91.187 GeV, 
nib = 4.5 GeV, sin 2 9w = 0.2315, and a = 1/128. We adopt a simple one-loop formula for 
the running strong coupling constant a s 

i \ oc a {m z ) . . 

a s(») = — 13=2^ 7 - v (3-1) 

where a s (mz) = 0.117 and rif is the number of active flavors at energy scale /i. In our 
numerical calculation to get the low energy scenario from the mSUGRA, the renormalization 
group equations(RGE's) [|T7| 1 are run from the weak scale Mz up to the GUT scale, taking all 
threshold into account. We use two loop RGE's only for the gauge couplings and the one- 
loop RGE's for the other supersymmetric parameters. The GUT scale boundary conditions 
are imposed and the RGE's are run back to Mz, again taking threshold into account. The 
decay widths of the intermediate Higgs bosons are considered at the tree level and these 
formula can be found in ref . ]18 ] . 



The neutralino pair production cross sections dependence on the factorization scale Q is 
illustrated in Fig.2, where we choose the mSUGRA parameters as m = lOOGeV, mi/2 = 
150GeV, A = ZOOGeV, tan/5 = 4 and /i > 0. The masses of x°i xl are hlAGeV and 98GeV '. 
For the X1X2 P arr production, the cross section changes from 28/6 to 30/6 when the scale Q 
changes from 0.2m £ o to m^o. For the X2X2 P arr production, the cross section changes from 
14/6 to 15.5/6 when the Q changes from 0.2m £ o to m £ o. When Q is larger than m £ o, the 
cross sections of these two processes are nearly independent of the factorization scale. 

In Fig. 3 we show the cross sections of gaugino-like neutralino pair production and the 
higgsino-like neutralino pair production. For the gaugino-like neutralino pair production, we 
choose the same values of mSUGRA parameters adopted in Fig.2. But in the higgsino-like 
neutralino case, we choose Mi = IbOGeV ', M 2 = 210GeV and /i = 90GeV in order to keep 
the masses of Xi an d X2 being the same as those in the previous gaugino-like neutralino case, 
the other parameters of the higgsino-like neutralino case at the weak scale are chosen being 
the same as those in the gaugino-like neutralino case. Then the difference of the cross sections 
between gaugino-like neutralino pair production and higgsino-like neutralino pair production 
only comes from the change of the coupling of x°i ~ Q — 1 and Xi ~X°j ~ H°(A°, h°, G°) where 
the matrix elements of N are changed. For the X1X2 P & ir production, the cross section in 
the higgsino-like case is smaller than that in the gaugino-like case, but for the X2X2 P a ^ r 
production, the cross section in the higgsino-like case is larger than that in the gaugino-like 
case. 

The cross sections for X1X2 an d X2X2 productions at the LHC versus the mass of X2 * s 
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shown in Fig. 4. We assume that Xi is the LSP and escapes detection. The input parameters 
are chosen as m = lOOGeV, A Q = 300GeV, tan (3 = 4, /i > and varying from 130GeV 
to 325GeV. In the framework of the mSUGRA, the masses of Xi an d X2 increase from AlGeV 
to 129GeV and from 81GeV to 251GeV respectively, and the masses of Higgs boson H° and 
A increase from 230GeV to 509GeV and from 22QGeV to 507GeV respectively with our 
input parameters. In our case, the mass of the Higgs boson h° is always smaller than 
m^o + m^o and the contribution to the cross section of h° is very small due to the s-channel 
suppression. Since the masses of H° and A are larger than m^o + m^o, the cross section 
will be strongly enhanced by the s-channel resonance effects when the total energy \/§ of 
the subprocess approaches the mass of H° or A . Since m^o > m^o, the pair production of 
X2X2 compared to X1X2 i s kinematically suppressed. For the X1X2 P & i r production, the cross 
section can reach 43 femto barn when = 130GeV, m^o = 41GeV and m^o = 81GeV. 
The cross section of X2X2 P a i r production can reach 23 femto barn when mm = 130 GeV. 
From the results of 0, we know that the cross sections of X1X2 anc ^ X2X2 P a i r production 
via quark-antiquark annihilation at the LHC are about 70 femto barn and 300 femto barn 
respectively, when the input parameters are chosen as m = 150 GeV, A = 300 GeV, 
tan/5 = 4, > and m^o = 81 GeV. The pair production rate of X1X2 Y ^ gluon-gluon 
fusion is about few ten percent of that via quark-antiquark annihilation at the LHC with 
the same input parameters. The cross section of X2X2 y i a gluon-gluon fusion is about few 
percent of that via quark-antiquark annihilation at the LHC. The neutralino pair production 
via gluon-gluon fusion at the LHC is of the same order of the next leading order(NLO) 
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QCD correction to the quark-antiquark production process. The neutralino pair production 
rate via gluon-gluon fusion should be considered in studying the pair pair productions of 
neutralinos at proton-proton colliders. 

In Fig. 5 we present the cross sections of neutralino pair productions versus tan/3 where 
the other four input parameters are chosen as m = 100 GeV, A = 300 GeV, /i > and 
m-i/2 = 150 GeV . When tan (3 increases from 4 to 32, the mass of Xi increases from 51 GeV to 
57 GeV and the mass of \\ increases from 98 GeV to 102.5 GeV . So the masses of neutralinos 
are nearly independent of tan f3. The masses of Higgs boson H° and A depend on tan f3 and 
decrease from 258 GeV to 165 GeV and from 254 GeV to 165 GeV respectively. The mass 
of h° is a function of tan (3 too, but keeps m h o < m^o + m^o. Since the couplings of Higgs 
bosons to quarks and squarks pair are related to the ratio of the vacuum expectation values, 
tan (3 should effect the cross sections substantially. The cross sections of the neutralino pair 
productions decrease when tan/5 goes from 4 to 6, and increase when 20 > tan/3 > 6 due to 
the effect of the coupling strength of Higgs bosons to quarks and squarks pair. For the X1X2 
pair production, when tanfi < 27, the masses of H° and A are larger than m^o + m^o, the 
s-channel resonance effects of H° and A on the cross section can be enhanced, while when 
tanf3 > 27, the masses of H° and A become smaller than m^o + m^o , then there will be no 
s-channel resonance effects and the cross section goes down rapidly. For the similar reason, 
the cross section of the X2X2 P a ^ T production increase in the range of 20 > tan (3 > 6 due 
to the effect of the Yukawa couplings of Higgs bosons to quarks and squarks pair and the 
s-channel resonance effects of the intermediate Higgs bosons, but it goes down again when 
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tan (3 > 20 where the masses of H° and A become smaller than m^o + m^o and the s-channel 
resonance effects of the Higgs bosons disappear. From the figure we can see that if we choose 
suitable input parameters, the cross section can be largely enhanced. Especially the cross 
section of the X1X2 P a i r production can even reach about 80 femto barn when tan {3 ~ 27. 

In our calculation we find that the discrepancy between our results from the newer sets 
of gluon densities in the MRST(mode 2) scheme and the older sets in the MRS(set G) 
scheme |TP[] is very small and cannot be distinguished in figures. 



4. Summary 

In this paper, we studied the pair production process of the neutralino via gluon-gluon 
fusion at the LHC collider. The numerical analysis of their production rates is carried out 
in the mSUGRA scenario with some typical parameter sets. The results show that the cross 
section of the neutralino pair-production via gluon-gluon fusion can reach about 80 femto 
barn for X1X2 P arr production and 23 femto barn for X2X2 P a i r production at the LHC 
collider. In some c.m.s energy regions of incoming gluons, where the s-channel resonance 
conditions are satisfied in the parameter space, we can see observable enhancement effects 
of the cross sections. We conclude that neutralino pair production via gluon-gluon fusion 
can be competitive with the quark-antiquark annihilation production process at the LHC 
and should be considered as a part of the NLO QCD correction to the quark-antiquark 
annihilation production process. 
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Figure Captions 

Fig.l The Feynman diagrams of the subprocess gg — > XiXj- 

Fig. 2 Dependence of the cross sections for the productions of X1X2 an d X2X2 pairs on 
the factorization scale Q with the mSUGRA parameters m = 100 GeV, mi/ 2 = 150 GeV, 
A = 300 GeV, tan/5 = 4 and /1 > 0. The masses of x°i x\ are 51.4 GeV and 98 GeV. 

Fig. 3 The cross sections for the neutralino pair production at the gaugino-like neutralino 
case and higgsino-like neutralino case at the LHC. For the gaugino-like neutralino case , 
the parameters are chosen in mSUGRA scenario with m = 100 GeV, mi/ 2 = 150 GeV, 
A = 300 GeV, tan/9 = 4 and /i > 0. For the higgsino-like neutralino case, we choose 
Mi = 150 GeV , M 2 = 210 GeV, /i = 90 GeV and the other parameters the same as those 
in the gaugino-like neutralino case. 

Fig. 4 Total cross sections of the process pp — > gg + X — > Xi Xj + A as function of m^o at 
the LHC. The solid curve is for X1X2 P arr production and the dashed curve is for X2X2 P & i r 
production. 
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Fig. 5 Total cross sections of the process pp — > gg + X — * Xi'X.'j + X as function of tan/5 
at the LHC. The solid curve is for X1X2 P 8 * 1 production and the dashed curve is for X2X2 
pair production. 
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Fig.(2) 
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7s = 14 TeV 



a(pp->gg + X-> X 2 X° 2 + X)[fb] 
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